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The irradiation of atoms by a strong IR laser field of frequency uj results in the emission of 
odd-harmonics of u ("IR harmonics") up to some maximal cut-off frequency. The addition of an 
XUV field of frequency qui larger than the IR cut-off frequency to the IR driver field leads to the 
appearance of new higher-order harmonics ("XUV harmonics") q ± 2K, 2q ± (2K — 1), 3q ± 2K, ... 
(K integer) which were absent in the spectra in the presence of the IR field alone. The mechanism 
responsible for the appearance of the XUV harmonics is analyzed analytically using a generalization 
of the semiclassical re-collision (three-step) model of high harmonic generation. It is shown that 
the emitted HHG radiation field can be written as a serie of terms, with the HHG field obtained 
from the three-step model in its most familiar context [P. B. Corkum, Phys. Rev. Lett. 71, 1994 
(1993)] resulting from the zeroth-order term. The origin of the higher-order terms is shown to be 
the ac-Stark oscillations of the remaining ground electronic state which are induced by the XUV 
field. These terms are responsible for the appearance of the new XUV harmonics in the HGS. The 
XUV harmonics are formed by the same electron trajectories which form the IR harmonics and 
have the same emission times, but a much lower intensity than the IR harmonics, due to the small 
quiver amplitude of the ac-Stark oscillation. Nevertheless, this mechanism allows the extension of 
the cut-off in the HGS without the necessity of increasing the IR field intensity, as is verified by 
numerical time-dependent Schrodinger equation simulation of a Xe atom shined by a combination 
of IR and XUV field. 

PACS numbers: 

03.65.-w, 42.50.Hz, 42.65.-Ky, 32.80.Rm 

Focusing intense linearly-polarized monochromatic infra-red (IR) laser pulses into gas of atoms can lead to the 
emission of high-energy photons with frequencies extending into the extreme ultraviolet (XUV) and X-ray region by 
high harmonic generation (HHG). All major features of HHG, such as its comb-like spectrum of odd-integer harmonics 
(to be called "IR harmonics"), its photons' maximal energy (to be called "IR cut-off') and the emission times of each 
harmonic, could be well reproduced using a semiclassical three-step (recollision) model 1, 2, 3J: under the influence of 
the intense laser field the electron of an atom tunnels out of the modified Coulomb potential, gains kinetic energy as a 
free particle in the field and finally may recombine with the parent ion to release the sum of its kinetic energy and the 
ionization potential as a high energy photon. The emission times of different harmonics are perfectly synchronized 
with the driver field, making the HHG process a promising method for the production of an adjustable coherent 
X-ray source. The current method of achieving the state of the art IR cut-off positions in the harmonic generation 
spectra (HGS) makes use of high-intensity few-femtosecond IR laser pulses. The main drawback of this method is that 
the electronic plasma which is inevitably formed at such high intensities, causes large dispersion on the propagating 
harmonics and severely limits their phase matching. The method to achieve higher-energy harmonics which will be 
presented here, doesn't suffer from this limitation, since it allows the usage of an IR source of moderate intensity 
which produces a small amount of plasma. It uses an XUV driving field which is shined on the atom simultaneously 
with the IR one. Taking the frequency of the XUV field larger than the IR-cut off frequency and inside some spectral 
window of the HHG generating gas, new higher order harmonics (to be called "XUV harmonics"), which were absent 
in the spectra in the presence of the IR field alone, could be produced, with frequencies well above the IR cut-off 
frequency, while the XUV photoionization could be suppressed, thus un-altcring the amount of electronic plasma. As 
will be shown, the main drawback of this method is that the new XUV harmonics have a relatively low intensity. 

The idea of contaminating the strong IR field with a second or more higher-frequency [usually ultra-violet (UV)] 
fields, is not a new one in the context of strong laser-matter interactions. The effect on the dynamical behavior of 
the electrons is dramatic, and such two-color (bichromatic) schemes have drawn a lot of attention in recent years. 
The additional field, if adjusted correctly, can induce stimulated emission [3], single-photon ionization J5j, or multi- 
photon ionization [6]. This ionization effect is actually utilized for attosecond (as) pulse-duration measurement, by 
the analysis of photoelectrons which are emitted from atoms exposed simultaneously to the as-pulses and strong IR 
field [TJ. The HGS obtained using bichromatic laser fields EO EH EH E2] , polychromatic fields E3E1ES], or even 
as-pulses [5] instead of the conventional monochromatic field, had been studied extensively as well. To the best of our 
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knowledge, in all above-mentioned studies, the frequency of the UV field was in the plateau of the HGS generated by 
the IR field alone (to be called here " IR HGS" ) , or close to its cut-off [16j EZ] • Such a UV field could only increase 
the efficiency of the existing IR harmonics and/or create new peaks (hyper raman lines) still in the support of the IR 
HGS, or maybe extend the IR cut-off to some limited extent [As will be shown here, it is the taking of the XUV field's 
frequency beyond the IR-cut off that pushes the cut-off position substantially to higher frequencies] . On the basis of 
the three-step (re-collision) model, it had been argued that the role of the UV field is to switch the initial step in the 
generation of high harmonics from tunnel ionization to the more efficient single UV-photon ionization |18j , or to assist 
the tunneling by transferring population to an excited state, from which the tunneling rate is larger |17[ll9j. especially 
if the UV photon energy matches some level transition EI] • This might explain the improved macroscopic HHG 
signal obtained in experiments: the UV-assisted ionization increases the number of atoms which participate in the 
HHG process and improves phase matching (the possibility that the improvement in HHG efficiencies is due to the 
interaction of the strong IR field with the created ions was shown to be implausible [22]). The above explanation 
doesn't apply, however, for the case that will be discussed in this paper, i.e. a case in which the high-frequency field in 
the bichromatic HHG scheme is in the XUV regime (and not in the UV one), with frequency well above the IR cut-off 
frequency. By choosing the frequency of the XUV field to fall inside some spectral window of the HHG generating 
gas, the XUV photoionization process could be eliminated (the same effect is achieved automatically for high-enough 
energies of the XUV photon since the single XUV photon ionization cross section scales as the 7/2-th power of the 
XUV photon's wavelength [23 ). Thus, in this case, the contribution of the XUV field to the generation of the XUV 
harmonics is not via affecting the ionization stage anymore. Another suggestion for the role that the high-frequency 
photons play in the HHG process was that they control the timing of ionization, and preferentially select certain 
quantum paths of the electron [23] . While this effect may lead to the enhancement of the low-order harmonics in the 
plateau, it can't account for the large enhancement in the cutoff and beyond (which is noticeable in the results in 
FigQ. 

A three-step model classical analysis of HHG suggests that the contribution of the XUV field to the kinetic energy of 
the returning electron is negligible. To see this, suppose we irradiate the atom with a linearly-polarized IR fundamental 
field of frequency ui, amplitude e™ and polarization (Ei(t) = e^e^ 1 cos(wi)) and an XUV field of frequency qoj 
(where q is a large enough number) and amplitude e™ (Eq(£) = ek£™ cos(gwi)) with the same polarization. By 
integrating the classical equation of motion while assuming that the electron is freed at time U with zero momentum, 
the following expression for the momentum of the electron is obtained: 

p(t)=Px(«)+Pq(t) (1) 

where Pi(i) is the momentum due to the IR field alone Pi(t) = ^j- [sin(u>t) — sin(ujti)] and Pq(i) is the momentum due 

to the XUV field alone p q (i) = ^j-[sin(qujt) — sin(qujti)] (e and m are the electron's charge and mass, respectively). 
The kinetic energy with two fields simultaneously present is: 

Ek{t) = ^ = EkAt) + E k A (t) + (2) 

where E k<1 (t) = ^ and E k) g(t) = Note that E k ^(t) oc (^) 2 and that the cross-term ^ e J_ are mucn 

smaller than E kt i(t) if e j^- < e™ . Usually, a large enough value of q will fulfill this condition, even if the IR and XUV 
fields have similar intensities. Thus, the additional XUV field will not affect the electron trajectories and will not 
contribute to their kinetic energy. For this reason the relative phase between the two fields doesn't play a role in the 
HGS, which is indeed verified in both classical analysis and quantum mechanical simulations (a small q, however, will 
affect the dynamics differently [§., I25|). In addition, assigning the electron a non-zero initial momentum to account 
for the photoelectric effect, will not increase its kinetic energy upon recombination. To conclude, the extension of the 
harmonic cutoff energy due to the inclusion of the XUV field (Fig(T]), isn't a result of an increase of the electron's 
kinetic energy upon recombination. 

Since the XUV field doesn't affect the kinetic energy of the electron trajectories (second step in the re-collision 
model), nor does it modify the ionization step (first step in the re-collision model), it must then influence the recom- 
bination step (third step in the re-collision model) of HHG. It is the purpose of this paper to prove this hypothesis. 
As will be shown later, the XUV field induces periodic ac-Stark modulations to the remaining ground electronic state, 
with the same frequency as the XUV field. The returning electronic wavepacket recombines with this modulated 
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ground state to emit the new XUV harmonics. The IR-HGS enhancement and the cut-off extension are a result of a 
single atom phenomenon, and not a macroscopic one. 

The paper is organized as follows: in section II we give the numerical results of the HGS obtained using a model 
Hamiltonian which describes a one-dimensional Xe atom subjected to a sine-square pulse of bichromatic field of 
frequencies uj and qui, for different values of q. In section III we briefly describe the semiclassical re-collision model 
of HHG, as is usually applied to the monochromatic case. We then modify the re-collision model to account for 
possible non-trivial time-dependence of the ground electronic state, and show that this modified re-collision model 
successfully reproduces the results presented in section II. In section IV we suggest an experiment based on the effect 
we discovered and conclude. 



XE ATOM DRIVEN BY A TWO-COLOR (ui, quS) LASER FIELD 



As an illustrative numerical demonstration of the IR cut-off extension in the HGS due to the addition of an XUV 
field, we studied a single electron ID Xe atom irradiated by a sine-square pulse supporting N oscillations of linearly 
polarized light of bichromatic field, composed of an IR laser field of frequency lu and amplitude e™ an d a high- 
frequency field of frequency qui and amplitude e™. The following time-dependent Schrodinger equation (TDSE) was 
integrated using the split operator method: 



ih^(x, t ) = S.^ + V (x)- exsm 2 [e? cos(u;t) + ef cos(quit)} W, t) (3) 

between the times < t < NT (T — — ) with atomic units {% = m = — e = 1) and with the wave function taken 
initially as the ground state <fci(x) of the field- free model Hamiltonian of a ID Xe atom, with the field-free effective 
potential Vq(x) = —0.63exp(— 0.1424a; 2 ). This potential supports three bound states, of which the two lowest ones 
mimic the two lowest electronic states of Xe, with energies —I p = — 0.4451a. u. and Ei — — 0.1400a. u.. The parameters 
used for the simulation were N = 50, ui = 0. 05695a. u. (A = 800nm), e" 1 — 0.035a. it. (corresponding to intensity of 
7|™ ~ 4.299 • 10 13 ^j), ef = O.OOOla.u. {If ~ 3.509 • W s W/cm 2 ). It should be noted that in this type of simulation 
the Born approximation is assumed: instead of solving two coupled differential equations, one for the evolution of the 
electron and the other for the propagation of the electromagnetic field, the electromagnetic field is assumed to remain 
unchanged during the interaction with the electron, and only a Schrodinger equation of motion for the electron is 
solved. 

In order to calculate the HGS the Larmor approximation |26j was assumed, and the time-dependent acceleration 
expectation value 



a(t) = ~ ^W>*)> + ^^ 2 (^F)[4"co S M) + 4 n cos(^)] (4) 

which is linearly proportional to the emitted field, was analyzed. The power spectra (HGS) of emitted radiation by 
the oscillating electron is proportional to the modulus-square of the Fourier-transformed time-dependent acceleration 
expectation value: 

f7(n) = |!| a( 0)|2 (5) 
where the acceleration in frequency space is given by the Fourier transform 



= m j Q a(t)e~ mt dt (6) 

Fignl shows the HGS for different values of q (q = 11, q = 25, q — 37^ and q = 52). The HGS in the presence 
of the IR field alone is also shown for comparison. Several features can be seen in the figure: The position of the IR 
cut-off (the cut-off of the HGS in the presence of the IR field only) is at the 15th harmonic where the maximal IR 
harmonic is the 29th one. When the high frequency field has a frequency still within the IR-HGS, its only influence 
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on the HGS is to modify the IR-harmonics (q = 11 in FigjT]). When the high frequency field has a frequency close to 
the IR cut-off, the IR-harmonics are modified, and new harmonics (above the 29th harmonic), which were not present 
with the IR field alone, appear (q = 25 in Fig(lJ which corresponds to an XUV radiation of wavelength 16nm). These 
new harmonics, that appear due to the addition of the XUV field only, will be termed XUV-harmonics. The position 
of the maximal XUV harmonic increases as the XUV frequency increases (q = 37 J|), where it becomes apparent that 
the XUV harmonics appear around the frequency of the XUV field. The number of new XUV harmonics reaches a 
maximal one (approximately twice the number of IR harmonics) whenever the value of q is either not odd integer or 
is greater than approximately twice the number of the maximal IR harmonic (q = 52 in FigjT]). 

In general, upon the addition of the XUV field of frequency qui to the IR field of frequency ui the harmonics q±2K 
(K integer) are either modified (if they were already present in the IR-HGS) and/or appear as new XUV harmonics. 
Moreover, the HGS possesses certain symmetries: with respect to its center at harmonic q, the distribution of the 
XUV harmonics is symmetric (i.e., for q = 37^, fi(33^a;) ~ a(Al-^ui), etc.) and upon variation of q it shifts but 
remains almost invariant. The structure of the HGS of the XUV harmonics (will be called XUV-HGS from now on) 
consists, in principal, of two new plateau-like regions and two new cut-off-like regions. For example, for q — 52 in 
FigjT] the harmonics of order 38-48 and 56-66 have a "plateau" character (constant intensity), and the harmonics 
32-36 and 68-72 have a "cut-off" character (constant phase, will be shown in Fig|5j. The XUV harmonics are 10-orders 
of magnitude weaker than the IR harmonics. 




FIG. 1: (color online) HGS obtained from a ID model Hamiltonian of Xe atom (Eq|3| irradiated by a 50-oscillation sine-square 
pulse of bichromatic laser field composed of a 800nm IR laser field of intensity I[ n ~ 4.299 • 10 13 H7 cm 2 and a 800/ij-nm XUV 
field of intensity IJp ~ 3.509 • W 8 W/crn 2 for different values of q: q — 11 (solid red line), q — 25 (solid green line), q — 37^ 
(solid blue line) and q = 52 (solid cyan line). The IR-HGS is shown as the dotted black line where the position of the IR cutoff 
is at the 15th harmonic. The addition of the XUV field could have several effects on the HGS, depending on its frequency with 
respect to the position of the IR cut-off: it could either modify the IR-HGS or slightly extend the IR cut-off (q = 11), or add 
new XUV harmonics to the spectrum (q — 25, 37^, 52). The XUV harmonics appear at q ± 2K, regardless of the value of q. 
The XUV-HGS possesses the following symmetries: with respect to its center q, the distribution of the new XUV harmonics is 
symmetric (i.e., for q = 52, a(48uj) ~ a(56u), etc.) and upon variation of q it shifts but remains almost invariant. 

Fig{2]shows the HGS for a larger intensity of the high-frequency field, ef = 0.0035a.u. {if ~ 4.299 • 10 n H7cm 2 ), 
where all other parameters are kept the same. For sake of clarity, only two values of q are shown q = 37 ^ and q = 52. 
As before, the set of new XUV harmonics q ± 2K (K integer) appears around q. It is 7-orders of magnitude weaker 
than the IR harmonics. In addition, an additional set of new XUV harmonics 2q ± (2K — 1) (K integer) appears 
around 2q and possesses the same above-mentioned symmetries: it consists of two new plateau-like regions and two 
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cut-off- like regions. The intensity of the XUV harmonics around 2q are, however, 14-orders of magnitude weaker than 
the IR harmonics. A third set of XUV harmonics 3q±2K appears around 3q. It possesses the same above-mentioned 
symmetries and is 21-orders of magnitude weaker than the IR harmonics. 
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FIG. 2: (color online) Same as in Fignl but for J|" ~ 4.299 ■ 10 11 W/cm 2 . For q = 37^ (solid blue line) three sets of XUV 



harmonics are observed: one at q ± 2K", a second one at 2q ± {2K — 1), and a third one at 3q ± 2K, all differ very much in 
their intensities but have the same general structure and the same symmetries pointed out in Fig jl] and in the text. For q = 52 
(solid cyan line) only 2 sets of XUV harmonics are observed, together with traces of the third set. The constant intensity of 
the IR-HGS (dotted black line) above the IR cut-off should be disregarded as it is due to numerical error and could be made 
as low as one wishes. 



For a larger intensity of the high frequency field this same trend is continued, giving rise to new sets of XUV 
harmonics around Aq ± (2K — 1), 5q ± 2K , etc. which could become nested but are well distinguished by the large 
differences in their intensities. Fig(3] and Fig|4j which show a top-view plot of the bichromatic HGS as function of 
q, suggest that the new sets of XUV harmonics emerge from the single set of IR harmonics. The sets of straight 
lines with different slopes in those plots correspond to different values of the integers n\ and n§ in the selection rules 
applicable for our bichromatic (cj, qui) scheme |27j : the possible harmonics that could be emitted are Q/uj = n\ + n^q 
(where n\ + rig = 2K — 1), or alternatively 
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(7) 



We shall symbolize the sets of lines according to their values of n\ and rig as (ni, rig). The set of lines parallel to the 
vertical axis of Fig [3] (infinite slope) correspond to rig — and odd values of m, i.e. to the IR harmonics f2 = (2K— 1)uj 
(in a perturbative picture, this corresponds to no absorption of qui photons, only an odd number of ui photons) and 
will be symbolized as (ni,riq) — (2K — 1,0). It is obvious that without the XUV field, only this set would appear in 
this type of figure. The next set of lines, with slope equal to unity, correspond to (m, rig) = (2K, 1), i.e. to the XUV 
harmonics ft — (q ± 2K)u that are shown also in FigJT] It is seen that these lines are much weaker than the lines 
corresponding to the set (m, rig) = (2K — 1,0). The strongest line in the set rig — 1 corresponds to (m, rig) — (0, 1), 
i.e. absorption of one quj photon and no u> photons. This line, together with the lines to its right [which correspond 
to (ni, rig) — (2, 1), (4, 1), ...] emerge from the same points on the function q = 1 form where the IR harmonic lines 
emerge. The lines (rii, rig) = (—2, 1), (—4, 1), ... are a mirror image of the lines (ni, rig) — (2, 1), (4, 1), with respect 
to the central line (n^n^) = (0,1). This hints us that the symmetry properties of the XUV-HGS (i.e., the fact 
that the distribution of the XUV harmonics is symmetric with respect to q, shifts but remains almost invariant upon 
variation of q, and consists of two "plateau" and two "cut-off" regions) are the result of the fact that the XUV 
harmonics are "born" from the same electronic trajectories that produce the IR harmonics. When the intensity of 
the high-frequency field is increased, as in FigjlJ also the sets (ni,riq) = (2K — 1,2) and (ni,ng) — (2K,3) appear, 
and they are 20-orders and 30-orders of magnitude smaller than the IR harmonics, as shown in Fig(2j Also these 
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XUV harmonics are "born" from the electronic trajectories that produce the IR harmonics. What is the physical 
mechanism leading to the formation of these new sets of XUV harmonics? 




FIG. 3: (color online) Top view of a log plot showing the HGS (red color-high intensity, blue color-low intensity) obtained from 
a ID model Hamiltonian of Xe atom (Eqj3| irradiated by a 50-oscillation sine-square pulse of bichromatic laser field composed 
of a 800nm IR laser field of intensity l[ n ~ 4.299 • 10 13 W/cm 2 and a 800/gnm XUV field of intensity if ~ 3.509 • VfW/ cm 2 for 
different values of q. The HGS shown in Fig[T]is obtained by taking cuts of the HGS shown here along the specific values of q. 
The lines in the spectrum are merely a manifestation of the selection rules given in Eq[7] Groups of harmonics corresponding 
to absorption of more than one XUV photon are absent, due to the weak intensity of the XUV field. That is, the only set of 
XUV harmonics that appear is the set q ± 2K. 




THE RE-COLLISION DESCRIPTION OF BICHROMATIC (u!,qui) HHG 



The process of HGS could be successfully described in terms of a very simple and intuitive model: the semiclassical 
re-collision model pQ . Let us describe, for the beginning, the process of HHG driven by an IR filed only. According 
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to this model, the electronic wavefunction at the event of recombination ^(r^ w t r ) could be described as a sum of 
the following continuum and bound parts 



*(r,tr) = il) b {r,t r ) + il) c {v,t r ) (8) 

It is assumed that the strong IR field ionizes the electron by tunneling from the initial ground state of the field-free 
Hamiltonian </>i(r), which is only slightly depleted during this process. It is assumed that the electronic wavefunction 
which remain bound, evolves under the field-free Hamiltonian only, i.e. accumulates a trivial phase only: 



#,() = f(r)e^ f (9) 

where —I p is the energy of the ground state. Under the strong field approximation, the freed electronic continuum 
part evolves under the external field only. Taking the direction of linear polarization as the x-direction from now 
on for simplicity, and assuming separability of the continuum wavefunction tj) c (r,t ~ t r ) in the x-coordinate and the 
2 other lateral coordinates for simplicity, the continuum wavefunction can be written as 



^ c (r,*w*r) = rp\{x,t™t r )^{y,z,tttt r ) (10) 
where the returning continuum part in the direction of polarization is some superposition of plane waves 

1 f 00 II E 

4(x,t r )=— dkip c \k,t r )e^ kx -^ (11) 
v2ir j-oo 

2 2 ~ | 

where k = fce x (k = |k|) is the momentum of the electron, = is the usual dispersion relation and ip c (k, t r ) 
are expansion coefficients which weakly depend on time. 

Using the total wavefunction at the event of recombination \I/(r, t « t r ), the time-dependent acceleration expectation 
value could be calculated. Keeping only the part which is responsible for the emission of radiation at frequencies other 
than the incident frequency w, the acceleration reads a(t) = ^ (^(r, t) | — VVb(r) |^(r, t)) r , where Vo(r) is the field-free 
potential. Assuming low depletion rate of the ground state (and hence, small population of the continuum wavepacket), 
the dominant terms that are responsible for the emission of radiation at frequencies other than the incident frequency 
u! are the bound-continuum terms 



a(t r ) = -$t(i(j b (r,t r )\ - VVb(r)|V'c(r,tr))r 
m 



(12) 



(the bound-bound term (^,(r,f)| — VVb(r)|-0b(r, i)) r is time-independent and doesn't radiate and the contribution 
of the continuum-continuum term is negligible) . After plugging the expressions in Eq 9JTT into Eq |12| it could be 
realized that the acceleration is composed of oscillating terms of the form 



~ (0) 2 2 

a(tr) = K/ dkij) (k,t r )e-^ I * +! ^t^ (13) 



■3C- 



where 



jj>° (k,t r ) = ---±=j c \k,t r ) I dMx^)VV Q {v)^{y,z,t r )e ikx (14! 



m V27T 

The emitted HHG field in a single re-collision event at t w t r is a burst of light which corresponds to the spectral 

continuum I p < hfl < I p + 3.17U P where 3.17C/ p = 3.17 4m ^ 2 is the value of the most energetic returning electron 
trajectory. In addition, since we assumed that the IR field pulls the electron along the x-direction, it is reasonable to 
assume symmetric evolution of the continuum wavefunction in the lateral plane, i.e., (y , z , t r ) — ip^(—y,—z,t r ). 



Since for atoms Vo(r) and <fii(r) are symmetric functions (and VVo(r) is antisymmetric), we get from Eq 14 that the 
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coefficient ip (k, t r ) has a nonzero component along the x-direction only, i.e., the acceleration a(t r ) points along the 
x-direction, as it should. 

When each single re-collision event is repeated every half cycle of the IR field, integer odd harmonics 2K — 1 are 
obtained in the HGS. To see this we compare two consecutive re-collision events at times t r and t r + %. Suppose 

we assume that ipc(x,t r ) was born at some initial time ti from 4>\{y). Therefore, since ipc(x,t r + ^) was born ^ 

after ipc(x,t r ), at the time ti + it was born from 4>i{Y)e^ Ip ^ since the bound state from which the continuum 

state tunnels out, has accumulated this phase. In addition, the two continuum functions ipi(x, t r ) and %j)\{x,t r + ^) 
are released in opposite spatial directions, because the IR field changes direction in two subsequent tunneling times. 
Therefore, the symmetry relation between tpc(x, t r ) and ipl(x,t r + -j) is 

iPUx,t r ) = 4(-x,t r + ^)e-* 1 ^ (15) 
and using the property = = E-k, we get by simple change of variables 

of the integral that 

$}{k,t r ) =i>}{-k,t r + T/2)e-^ +E "^ (16) 
Using the relations mentioned before (Vo(r) = V (— r), </>i(r) = <^i(— r), VVo(r) = — VVo(~ r )i V^(2A-Mr) = 



tp^(y, z, t r + T/2) = t/j^(—y, —z, t r )) the following symmetry is obtained from Eq 14 

i> ( °\k,t r ) = ~^°\-k,t r +T/2)e-K^ +Efe) ? (17) 
Using this symmetry, we get from Eq |13| that: 



a(t r + T/2) = -a.(t r ) (18) 

The acceleration vector is periodic in T and alternates directions between subsequent re-collision events, i.e., its only 
nonzero components in its Fourier expansion correspond to odd integer harmonics of lu. This is the origin of the well- 
known odd selection rules of monochromatic HHG. We shall therefore symbolize the acceleration which is responsible 
for the emission of odd harmonics (Eq{T|]) as a(°)(i). 

The three-step model described above assumes that the only time evolution of the remaining bound part of the 
electronic wavefunction is to accumulate a trivial phase, as given in Eq(9] This assumption, however, describes only 
the leading term in the time evolution of the bound part. In reality, due to the ac-Stark effect induced by the IR field, 
the electron adiabatically follows the instantaneous ground state of the potential which periodically shakes back and 
forth by the IR field. If one carries out a TDSE simulation and looks at the electronic wavefunction in the field-free 
potential region during the action of the IR field, one sees that it oscillates back and forth with the same frequency 
of the IR field. The time evolution of the bound part ?/>b(r,i) should be therefore corrected from the trivial one 
given in Eq(9] In addition, for common field intensities, the ac-Stark correction to the instantaneous ground state 
energy is negligible, and we may therefore assume that the instantaneous ground-state energy is almost constant (I p ). 
More importantly, the ac-Stark effect induces a periodic motion of the wavefunction as a whole, without deforming 
it. Relying on these facts, we approximate the instantaneous ground state wavefunction as 



V> 6 (r, t) = 4>i{x + ef* cos(wt), y, z)e+^ I " t (19) 

We have assumed the simplest time-dependence in ipb(^, t) that would still give periodic modulations at frequency ui. 
It should be noted that in the language of non-Hermitian quantum mechanics [28 , 3CJ] this expression approximately 
describes the resonance Floquet state which evolves from the ground state 4>i(r) upon the switching of the IR field. 

The quiver amplitude e° ut of the spatial oscillations of the ground state is of the order of e° ut = ( £2+ j 1 - )2 _ a) 2 (this 
is approximately the quiver amplitude of a an electron bound in a short-range potential of the type used here, driven 
by an IR field of amplitude £™), i.e., a tiny fraction of a Bohr radius, provided that the laser's frequency doesn't 
match some level transition. The bound part may therefore be expanded in a Taylor serie as 
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Mr, t) = e+i^lM?) + 4 Ut cos(wt)^0i(r)} 



(20) 



Calculation of the time-dependent acceleration expectation value using the total wavefunction at the event of 
recombination ^(r,t rs t r ) with the modified bound part ipb(r,t r ) given in Eq 20 and keeping again only the terms 
that are responsible for the emission of radiation at frequencies other than the incident frequency w (the bound- 
continuum terms), yields: 



a(i r ) = aS°\t r ) + cos(ut r )a w (t r ) 



(21) 



where a^(t r ) is given in Eq 13 and 



rOO ~(1) 2 2 



(22) 



where 



5,(1) 

1p (Mr) 



1 



m J2-K 



dx 



VV (r)^(y,z,t r )e 



ikx 



(23) 



We see that the inclusion of the Stark effect contributes a new term e\ ut cos(wt r )a' 1 ^(t r ) to the acceleration. In a 
single re-collision event at t « t r , this term produces two bursts of light. One corresponds to the spectral continuum 
I p + hui < fril < I p + 3.17U P + fiuu and the other to I p — hu < ~h£l < I p + 3.17U P — Tiuj. These bursts of light are, 

however, much weaker than the one which results from a^(t r ), since, as we recall, the factor e° ut = (g, + j 1 )2_ [ ^2 is 
small for common IR laser frequency and intensity (in the context of HHG experiments). It can be said in general 
that each electron trajectory (plane wave) with kinetic energy Ek, recombines with the nucleus to emit radiation of 
energy I p + Ek, and two "duplicate" photons with energies I p + Ek — fioj and I p + Ek + Tilo, at the same emission times. 
In a multi recollision sequence, only some of these photons will appear in the HGS, as dictated by selection rules 
which we are about to prove. In addition, since the functions and ip^r(y, z,t r ) are symmetric with respect to 

that the coefficient tp (k,t r ) and the acceleration a 1 (f r ) point 



22 



2:-! 



y, z and VVo(r) = — VVb(— r), we get from Eq 
along the x-direction, as they should. 

When each single re-collision event is repeated every half cycle of the IR field, the new term in the acceleration 
contributes integer odd harmonics 2K — 1 to the HGS. To see this we compare two consecutive re-collision events at 
times t r and t r + ? and note that the symmetry relation between tpc(x,t r ) and il>c{x, t r + -j) from Eq 

symmetry relation between ip c \k,t r ) and tp c (fc, t r + from Eq 



1G 
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and the 



still hold since the modification of the bound 



part of the electronic wavefunction has no influence on the continuum part 

Using the facts tha 
obtained from Eq(23] 



Using the facts that d^ r ) and VVb( r ) are antisymmetric with respect to inversion of x, the following symmetry is 



Tjj (1 \k,t r ) = +<ip il (-k,t r + T/2)e-* {Ip+Ek) % 



(24) 



Using this symmetry, by calculating a^\t r + T/2) using the definition given in Eq 22 together with the relation 
given in Eq(24j we get that: 



(t r + T/2) = +a^(t r ) 



(25) 



The acceleration vector '(t r ) is periodic in T/2, i.e., its only nonzero components in its Fourier expansion 
correspond to even integer harmonics of 10. The term which is responsible for the emission, e° ut cos(uit r )a^ (i r ), 
however switches signs every T/2, therefore giving rise to odd harmonics in the HGS (cos a cos (3 — \[cos(a — (3) + 
cos(a + /?)])• That is, the radiation resulting from the ac-Stark oscillations of the bound electron is composed of odd 
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integer harmonics of w, like the radiation which results from a^°\ The two fields, emitted by a(°) and e° ut cos(wi T .)a( 1 \ 
interfere with each other in general. However, since the field resulting from a' 1 ) is much weaker, it is completely masked 
by the field produced from a( ). 

The effect of the ac-Stark oscillations on the HGS could be summarized as follows: In a single re-collision event 
each electron trajectory (plane wave) with kinetic energy Ek , recombines with the nucleus to emit radiation at energy 
Ip + Ek, and also, due to the ac-Stark effect, two weaker "duplicate" electromagnetic waves with energies I p + Ek — noj 
and I p + Ek + nuj, at the same emission time. In a multi re-collision sequence, due to the symmetry properties 
discussed above, only odd-harmonic photons will appear in the HGS. The contribution to each plateau odd harmonic 

in the HGS comes, in principle, from six emission times: the first corresponds to the recombination of the short 
trajectory with kinetic energy h£l — I p , the second corresponds to the "duplicate" recombination resulting from a 
different short trajectory, with kinetic energy h£l — huj — I p (the ac-Stark oscillations of the ground state at frequency 
uj will make the final energy of the emitted photon Ml — hui — I p + Tiuj = tiil — I p ) and the third corresponds to 
the "duplicate" recombination resulting from a different short trajectory, with kinetic energy hfl + Klu — I p . Three 
additional emission times result from three long trajectories in the same manner. Because of the large differences in 
intensities, usually only the two "usual" emission times attributed to the short and long trajectories at kinetic energy 
Ml — I p will contribute. 

The effect of the ac-Stark oscillations on the HGS in this example is not large, since the harmonics produced by 
this mechanism are completely masked. However, high enough frequency of the ac-Stark oscillations, well above the 
IR cut-off, will cause the appearance of high energy photons which were not present at all in the IR HGS. In this case, 
the lack of contribution of the ordinary re-collision mechanism (neglecting the ac-Stark effect) to the appearance of 
these new high-energy photons makes the ac-Stark effect the only important one. 

The way to induce ac-Stark oscillation of high frequency is by applying a second high-frequency XUV field, in 
addition to the IR one. The emphasis is that in order to see the new harmonics produces by the "duplicate" 
trajectories due to the ac-Stark effect, the second field should be close to the IR cut-off or above, otherwise the new 
harmonics will be masked by the already existing IR ones. 

Suppose then that we shine the atom with an IR field of frequency oj and an XUV field of frequency quj. As 
discussed before, the XUV field, provided that it has a large enough frequency, doesn't affect the electron trajectories. 
Its only influence is therefore on the recombination process. This field induces ac-Stark oscillation in the exact same 
way as the IR field did. We may assume the same approximations used before and approximate the instantaneous 
ground state wavefunction as 

ip b (r, t) ^<p l {x + el ut cos(wi) + ef l cos(qut),y, z)e + ^ I " t (26) 

The quiver amplitude e^ ut of the spatial oscillations of the ground state is even smaller than e^ u± because of the 

high frequency of the XUV field ( e? ut = (g 2+ /)2_^2^,2 — ^js- The bound part may be expanded in a Taylor serie 
as before 



cos(u>t) 



d 

; cos(qwi)] — <Ai(r)} 



(27) 



and the time-dependent acceleration expectation value is calculated using the total wavefunction at the event of 
recombination ^(r,t « t r ) with the modified bound part ipi,(r,t r ) given in Eq 27 Keeping again only the bound- 
continuum terms, we get: 



a(t r ) 



= a<°> 



(tr) + [ef* cos(Lut r ) + ef l cos(qcot r )]a^ (t r ) 



(28) 



where aS°\t r ) is given in Eq[l3knd a^(t r ) in Eq 
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The ac-Stark effect at the frequency of the XUV field contributes a new term e- cos(qu>t r )aS ' (t r ) to the acceler- 
ation. In a single re-collision event at t « t r , this additional term produces two bursts of light. One corresponds to 
the spectral continuum I p + qhcu < Ml < I p + 3.17U p + qhuj and the other to I p — qhu < Ml < I p + S.17U P — qhuj 
which, in case that the XUV field is well above the IR cut-off, could be written as qhuj — I p < Ml < qfouj — I p — 3.17U P . 
These bursts of light are much weaker than the one which results from am(t r ), nevertheless they have a significant 
impact on the HGS since they are the only source of new XUV harmonics which appear now in the HGS. As explained 
before, each electron trajectory (plane wave) with kinetic energy Ek, recombines with the nucleus to emit radiation 
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of energy I p + Ek, and at the same time two "duplicate" photons with energies I p + Ek — qhu — qhuj — I p — Ek and 
Ip + E k + qhu. 

When each single re-collision event is repeated every half cycle of the IR field, the new term in the acceleration 
contributes to the HGS integer even harmonics around q, i.e. q ± 2K. To se e t his we compare two consecutive 

aW(t r 
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still holds here: 



T/2) 



re-collision events at times t r and t r + ^ and note that the result obtained in Eq : 

+aW(i r ). The acceleration vector a^(t r ) is periodic in T/2, i.e., contributes even integer harmonics of lo. The term 
which is responsible for the emission, e^ ut cos(q(jjt r )a^' (t r ) therefore gives rise to the appearance of the harmonics 
il = (q±2K)uo in the HGS (cosacos/3 = |[cos(a— j3) +cos(a+/3)]). That is, the radiation resulting from the ac-Stark 
oscillations of the bound electron is composed of even integer harmonics of lo around the harmonic q. This field is 
much weaker than the harmonics produced by a^ ). In case that q is an odd integer, the term e~ ut cos(qoot r )a^ (t r ) 
will produce odd harmonics, which will be masked by the IR harmonics produced from a*- ) , provided that q is well 
below the IR cut-off harmonic. However, in case that q is close to or above the IR cut-off harmonic or is not an odd 
integer, new harmonics (defined before as "XUV harmonics"), which were not present in the HGS in the presence of 



the IR field alone, will appear. These harmonics will be 



) 2 -times weaker than the IR harmonics, and 



could therefore be distinguished from the IR harmonics by their intensity. Since these harmonics are originated by the 
same electronic trajectories which produce the IR harmonics, their emission times are correlated with the ones of the 
IR harmonics. In a single re-collision event each electron trajectory (plane wave) with kinetic energy Ek, recombincs 
with the nucleus to emit radiation at energy I p + Ek, and also, due to the ac-Stark effect, two weaker "duplicate" 
electromagnetic waves with energies I p + Ek — qhoo and I p + Ek + qhu, at the same emission time. This correlation is 
kept also in the multi re-collision process and is manifested in the HGS: the structure (amplitude and phase) of the 
XUV harmonics (XUV-HGS) = (q±2K)u> is derived from the structure of the IR-HGS. This is true for every value 
of q but could be most easily seen if q is well above the IR cut-off harmonic, since in this case the XUV harmonics 
are separated and are not nested in the IR-HGS. If we look at the case q — 52 in Fig[l] we see that the structure 
of the XUV-HGS between the orders 54 and 74 resembles the structure of the IR-HGS between the orders 1 and 21 
(54 — 52 — 1 = 1 , 74 — 52 — 1 = 21) . In addition, within the XUV-HGS, the structure of harmonics between the orders 
32 and 50 is a mirror-image (with respect to the 52-nd harmonics) of the structure of harmonics between the orders 
54 and 72. That is, the XUV-HGS consists of two new plateau-like regions (harmonics of order 38-50 and 54-66), 
derived from the same electronic trajectories which form the IR-HGS plateau (harmonics of order 1-15), and two new 
cut-off-like regions (harmonics of order 32-36 and 68-74) , derived from the same electronic trajectories which form the 
IR-HGS cut-off harmonics (harmonics of order 17-23). This structure of the XUV-HGS is invariant to the value of q. 

What happens if we now take higher-order terms in the Taylor serie expansion of the ground state wavefunction in 
Eq |26| ? Suppose we take also the second-order term in the Taylor serie and expand the bound state according to 



■0b(r,i) 



{Mr)- 



lef* cos(^) + e e»* cos(g^)]^0 1 (r) + [e°«* cos(u,t) + e$* cos(g^)] 2 ^0 1 (r)} (29) 



1 dx 



2dx 2 



The time-dependent acceleration expectation value, keeping again only the bound-continuum terms, would read: 



a(t r ) = a (0) (i r ) + [ef*coe(wtr) + ef l cos(^t r )]a (1) (t r ) + [e° ut cos(^ r ) + ef l cos(qujt r )} 2 a {2) (t r ) 
= a< >(tr) + [ef* cos{ujt r ) + e? ut cos(gcO] a W(t r ) 



+ <! ( £ r*) 2 



1 + cos(2o;t r .) 



~out out 

2 1 « 



cos[(<7 + l)ujt r ] + cos[(g — l)wt r 



aS°'{t r ) + [ef * cos(Lut r ) + £? ut cos(qujt r )}a^> (t r ) 
(er 4 ) 2 + (e°g ut ) 2 (ei ut ) 2 



1 + cos(2qu>t r ) 



a (2) (M 



c OUt c-OUt 

cos(2o;t r ) H — - — cos[(g + l)u}t r ] 1 



cos[(g — l)wt r ] H — cos(2qujt r ) \a^ 2 



where aS°\t r ) is given in Eq 13 aW(t r ) in Eq 22 and a( 2 )(i r ) is 



x(2) 



aW(t r )=K/ dk^ (k,t r )e-^ [I " + '^ ]t 



(31) 



where 
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f\k,t r ) = -~^i,}{k,t r ) [ X d 3 r d2 p^WVo(r)^(y,z,t r )e ikx (32) 

m V 27T J-oo C>X 



The new term +[e° u * cos(ujt r ) + e~ ut cos(qtjt r )] 2 a^ 2 - ) (< r ), which results from the inclusion of the second-order 
term in the Taylor serie expansion of the bound wavefunction, produces 10 weaker "duplicate" bursts of light 
in a single recollision event at t ~ t r . This is due to the fact that the term multiplying a( 2 )(i r ) has 5 differ- 
ent frequency components. These bursts of light are much weaker than the one which results from a^°\t r ) or 
+ [e° ut cos(wt r ) + cos(qujt r ))a^ (t r ) . Nevertheless, it will be shown immediately that the last two bursts (resulting 

from -Jj — cos(2qu)t r )a^ (t r )) have a significant impact on the HGS since they are the only source of new XUV 
harmonics which appear around the harmonic 2q. 

Let us analyze now what happens when each single re-collision event is repeated every half cycle of the IR field . It 

the 
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was shown in Eq 18 that the term aS°\t r ) contributes odd harmonics f2 = (2K — l)u> to the HGS. From Eq 
term e° ut cos(u)t r )a^ (t r ) also contributes odd harmonics = (2K — 1)uj (with relative amplitude of the electric field 
of e\ ut ) but the term eq ut cos(qu>t r )a^ (t r ) contributes the XUV harmonics Q = (q ± 2K)u. Under the assumption 

that </>i(r) is symmetric, 9 is also symmetric. Therefore, in complete analogy to what was shown in Eq 



18 



it 

can be shown that the term a^ 2 \t r ) contributes odd harmonics of to since in two consecutive re-collision events at 
T 

2 



times t r and t r + % the following symmetry holds 



a^iU + T/2) = -a^(t r ) (33) 

The term which is responsible for the emission, [e° ut cos(ujt r ) + e~ Mt cos(<7wi r )] 2 a( 2 )(t r ) therefore gives rise to the 
appearance of the following five sets of harmonics (with the amplitudes given in parentheses): f2 = [0 ± (2K — l)]ui — 
(2K-l)w [(sf ut ) 2 + (e° s ut ) 2 ],Q= [2±(2K-l)]w = (2K-l)w [(si ut ) 2 ], il = [q + l±{2K-l)]u = (q±2K)u [e^ef**], 
n = [q - 1 ± (2K - l)]w = (q ± 2K)u [ef^ef*], = [2q ± (2K - [(s° 9 ut ) 2 ]. Out of these five sets, the first 
four contribute harmonics which are masked, due to their low intensity, by IR harmonics or by the XUV harmonics 
resulting from the lower-order term £~ u * cos(ga;t r )a( 1 ^(i I .). The fifth set, however, contribute a new set of harmonics: 

Q = [2q± {2K — l)]w. This set is |(e^"*) 2 | 2 — (;pjja ) 4 -times weaker than the IR harmonics, or |e?"*| 2 -times weaker 
than the set f2 = (q ± 2K)li of XUV harmonics and could therefore be distinguished from these set of harmonics by 
their intensity. As before, these harmonics are originated by the same electronic trajectories which produce the IR 
harmonics, and their emission times are correlated with the ones of the IR harmonics. This is manifested in the HGS: 
the structure (amplitude and phase) of the new set of XUV harmonics f2 = [2q± (2K — l)]u> in the XUV-HGS is 
derived from the structure of the IR-HGS. It consists, like the set f2 = (q ± 2K)uj, of two new plateau-like regions, 
derived from the same electronic trajectories which form the IR-HGS plateau, and two new cut-off-like regions, derived 
from the same electronic trajectories which form the IR-HGS cut-off harmonics. This structure of the XUV-HGS is 
almost invariant to the value of q, as can be easily seen from Fig(2] 

By generalizing the procedure described before and taking higher-order terms in the Taylor serie expansion of 
the bound wavefunction, it is apparent that the n-th order term contributes a new set of XUV harmonics CI = 

[nq ± (2K — 1 + mod(n, 2))]lu which is |(e^ ut ) n | 2 — ( ^7U2- ) 2rl -times weaker than the IR harmonics. Note that indeed 



in Fig 

( 



for q — 52, the intensity of the plateau harmonics of the set = (q ± 2K)co is indeed l^"')! 2 — (^t) 2 



0.0035 \2 



q 1 1 yq-'uj 2 



52 2 0.05695 2 



) 2 ~ 1.6 • 10 7 -times weaker than the IR plateau harmonics, and the intensity of the plateau harmonics 
is indeed |e^ ut )| 4 — 2.5 • 10~ 14 -times weaker than the IR plateau harmonics, in 



of the set Cl = [2q ± (2K - 1) 

agreement with our theory. Fig 5j which plots the intensity of different harmonics as function of the amplitude £ ~" of 
the XUV driver field, shows indeed that the intensity of harmonics from the set £1 = (q ± 2K)oj scale quadratically 
with while harmonics from the set = [2q± (2K — l)]ui scale as (e™) 4 - The obtained sets of XUV harmonics 
ft = [nq ± (2K — 1 + mod(n, 2))]uo are exactly as predicted by the selection-rules given in EqjT] Hence, the inclusion 
of higher-order terms in the Taylor serie expansion of the bound wavefunction, leads to the generalization of the 
scmiclassical three-step model. This allows us to obtain the selection rules for the high harmonics which are obtained 
upon the addition of an XUV field to an IR one, which are in complete agreement with the ones obtained using 
Floquet theory [27]. Moreover, the intensities of the XUV harmonics can also be quantified. 

The above theory predicts that both the IR and the XUV harmonics are emitted by the same IR trajectories and are 
therefore correlated in their emission times. This could be easily verified by plotting the time-frequency distribution 



13 





17 
—11 

IX - 




—la 
-19- 




-20 - 


a 


-21 - 













_?-> - 








-23- 




-24- 




-25- 




-26- 



■O=74o) 

Q=97(0 









-2.3 -2.25 -2.2 -2.15 

FIG. 5: (color online) log-log plot of the intensity of different harmonics a(Sl) [Q = 60 (solid red line), SI — 74 (dashed 
green line), = 89 (dash-dotted blue line), SI — 97 (dotted cyan line)] obtained from the ID model Hamiltonian of Xe atom 
(Eq|3| irradiated by a 50-oscillation sine-square pulse of bichromatic laser field composed of a 800nm IR laser field of intensity 
I\ n ~ 4.299 • 10 13 W/cm 2 and a 800/52-nm XUV field, as function of the XUV field's amplitude ejg. This intensity varies 
between 0.005a.«. (ij!J ~ 8.773 • W^W/cm 2 ) and 0.007a.u. (j|2 ~ 1.719 • 10 12 W/cm 2 ). All 4 graphs are linear. According to 
the lower frame of Figj2] (q — 52), the harmonics SI = 60 and SI = 74 belong to the set q ± 2K of XUV harmonics, and as 



such should depend on £52 as o(Sl = q ± 2K) oc (£52) (see Eq 28 1. Indeed the slopes of the graphs corresponding to SI = 60 



and S2 = 74 are 2.000 and 2.0001, respectively. The harmonics fi = 89 and SI = 97 belong to the set 2q ± (2K - 1) of XUV 
harmonics, and as such should depend on £52 as a [SI = 2q ± (2K - 1)] oc (£^) 4 (see Eq|30|. Indeed the slopes of the graphs 



corresponding to SI = 89 and SI = 97 are 4.0045 and 4.0015, respectively. This check confirms the correctness of our theory in 
this range of parameters. 

of high harmonics, i.e., by analyzing the Gabor-transform (windowed Fourier transform) of the acceleration instead 
of the Fourier transform: 

-1 pNT 2 

G ^ n ' t o)^^fJ Q a{t)e- { ^e- lM dt (34) 

where r is the window's width. This analysis yields a mixed time-frequency signal, i.e., not only the frequency 
components appearing in the acceleration but also their time of appearance to. Fig(6] shows the time-frequency 
distribution of high harmonics (a (ft, to) = ^g\G a {Q, to)\ 2 ) obtained from the time-dependent acceleration expectation 
value whose spectra is given in Fig(2]for q — 52, for the times 26T < t < 21 ST . In accordance with the semiclassical re- 
collision model, different harmonics are emitted repeatedly every half cycle. In accordance with the theory developed 
here, the time-frequency distribution of the new sets of XUV harmonics matches the one of the IR harmonics, including 
the reflection symmetry of these sets (around q,2q...). The most visible demonstration of this property is shown for 
the IR cut-off harmonics (the 15th-29th harmonic), who are emitted at times t ~ 0.710T + 0.5nT, in accordance with 
the semiclassical re-collision model. At those instants, also the 32nd-36th and the 68th-74th harmonics, which are 
the XUV cut-off harmonics of the set q± 2K, and also the 85th-87th and the 121st-125th harmonics, which are the 
XUV cut-off harmonics of the set 2q± (2K — 1), are emitted. They are hence produced by the IR cut-off trajectories. 

SUGGESTED EXPERIMENT AND CONCLUSIONS 

The mechanism described here raises the question whether XUV harmonics should be self-produced in any 
monochromatic HHG experiment, as high-harmonic radiation generated by the leading edge of the IR pulse, co- 
propagates with the IR field to form a bichromatic driver field in the last part of the medium. The answer is that the 
high-harmonic radiation generated in the medium is too weak to considerably modify the IR-HGS. If, on the other 
hand, a stronger high-frequency source is used, and in particular with frequency higher than the IR cut-off frequency, 
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FIG. 6: (Color online) Top view [pink (dark gray) colorhigh intensity, ye llow 
square of the Gabor-tansformed acceleration expectation value 
described in Fig[2]for q = 52, as function of to and Q. 



Eq 34 



(bright gray) colorlow intensity] of the absolute 
= 0.1T) of the quantum mechanical simulation 



new XUV harmonics should appear. This source could be a free-electron laser, but could also be a HHG-based 
source. Generation of HHG pulses of the 27th harmonic of Tksapphire laser, with width of 30fs and output energy 
of 0.33/^J per pulse, has been shown to be feasible in Ar |29j . When such an harmonic pulse is focused to an area 
of (10^m) 2 , intensities of J^? — 10 13 W/cm 2 may be reached. It isn't unreasonable to assume that higher-harmonics 
couldn't be generated with similar output energies. Even a reduction of 5 orders of magnitude in the intensity of 
the high-harmonic field will make the effect still visible. For example, by generating the 45th harmonic in Xe or He, 
filtering it out of the HGS and focusing it into a jet of Kr together with an IR field that is sufficient to generate IR 
cut-off at the 19th harmonic or so, new high harmonics should appear. Kr has a spectral window between photon 
energies of 50eV (the 31st harmonic) and 90eV (the 57th harmonic). Therefore, shining a 45th-harmonic field on it 
will not cause singlc-XUV photon ionization and will cause the appearance of new XUV harmonics in this spectral 
window, without the necessity to increase the intensity of the IR field. 

Alternatively, in case that the frequency of the high-harmonic field isn't large enough, one can reduce the intensity 
of the IR field in order to decrease the IR cut-off. The key point here is that in order to see the effect of appearance 
of new XUV harmonics the XUV field should be strong enough and of frequency higher than the IR cut-off frequency. 
Above all, the gas in which the bichromatic HHG is generated, should be transparent in some spectral band around 
the frequency of the XUV field, otherwise the generated XUV harmonics and/or the seed XUV field will be absorbed. 

In conclusion, we have shown that the addition of an XUV field to a strong IR field leads to the appearance of new 
harmonics in the HGS. The results of the semiclassical analysis and the quantum numerical simulations suggest that 
this is a single-atom phenomena. The XUV field induces ac-Stark modulations on the ground state and affects the 
recombination process of all returning trajectories, and leads to the generation of higher harmonics whose emission 
times and intensities are well related to the ones of the harmonics in the presence of the IR field alone. Using this 
mechanism, harmonics with unprecedented high frequencies could be obtained in HHG experiments. According to 
our mechanism, the emitted HHG radiation field could be written as a serie of terms, with the zeroth-order term 
equal the HHG field which is obtained from the three-step model in its most familiar context pQ. The higher-order 
terms become important when, in addition to the IR field, an additional XUV field is shined on the atom since they 
solely are responsible for the appearance of new harmonics in the HGS. 
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